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ABSTRACT 

Aims. The role of HeH + has been newly assessed with the aid of newly calculated rates which use entirely ab initio methods, thereby 
allowing us to compute more accurately the relevant abundances within the global chemical network of the early universe. A compar- 
ison with the similar role of the ionic molecule LiH + is also presented. 

Methods. Quantum calculations have been carried out for the gas-phase reaction of HeH + with H atoms with our new in-house code, 
based on the negative imaginary potential method. Integral cross sections and reactive rate coefficients obtained under the general 
conditions of early universe chemistry are presented and discussed. 

Results. With the new reaction rate, the abundance of HeH + in the early universe is more than one order of magnitude larger than in 
previous studies. Our more accurate findings further buttress the possibility to detect cosmological signatures of HeH + . 
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1. Introduction 

The existence of HeH + in astrophysical environments has 
been extensively discussed in several papers. Roberge & 
Dalgarno (1982) surmised that HeH + is produced in planetary 
nebulae and in dense molecular clouds under the presence of 
X-Ray and extreme UV ionization sources. Liu et al. (1997) in- 
vestigated the possible detection of HeH + in NGC 7027, a very 
high electron density nebular object, via its j = 1-0 strong ro- 
tational line at 149.14 /mi, concluding that its identification is 
complicated by the accidental near-coincidence of that transition 
frequency with the 149.09 and 149.39 jum lines of CH. 

In the last years the presence of HeH + in low temperature 
and low density helium-rich white dwarfs has further been in- 
vestigated: Harris et al. (2004) demonstrated that this molecule 
is the dominant positive ion in such objects and strongly affects 
their opacity. Engel et al. (2005), finally, showed that it could 
also be possible to detect HeH + in cool helium stars as those 
studied by Saio et al. (2000). 

Additionally, the discovery of the very metal-poor stars 
(Christlieb et al. 2002, Frebel et al. 2005) strengthened the quest 
for Population III stars. It then follows that to study the evolu- 
tion of primordial molecules formed during post-recombination 
era is a fundamental step to gain a better understanding of the 
"how and when" such stars were formed. 

To establish the role of the main molecular species (LiH, 
LiH + , HeH + , H 2 , HD + , HD) during the gravitational collapse 
which led to the first stars formation, it is however well-known 
that their opacity data are indeed needed. It is also just as impor- 
tant to know the elementary mechanisms which lead to the for- 
mation and destruction of such molecules and the corresponding 



rates that were present at the conditions of the astrophysical en- 
vironment. In other words, we have to know as much as possible 
the chemical history of all the above species. 

An extensive analysis of the lithium chemistry has been re- 
ported in a recent paper (Bovino et al. 2011) while an analo- 
gously accurate study of the reactions which involve helium- 
based molecular species, and in particular the HeH + molecular 
cation in reaction with H, is still missing. The latter is indeed 
considered to be the first molecule to appear in the Universe 
(Lepp et al. 2002): in an environment dominated by a low- 
density regime where 3 -body collisions are unlikely to occur, the 
radiative association is thus the main path to form that molecule. 
Hydrogen and helium, which are the most abundant species 
which would have existed at the early universe epoch, can there- 
fore combine to form HeH + by the familiar reaction: 



He + H + 



HeH + + hv. 



(1) 



for which the rate coefficient has been obtained from Zygelman 
et al. (1998) and fitted by Stancil et al. (1998). On the other hand, 
the main destruction path for the formed HeH + molecules is still 
the reaction with H: 



HeH + + H 



He + H£, 



(2) 
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which however has never been studied with sufficient accu- 
racy from a theoretical point of view. Crossed-beam experimen- 
tal data for energies larger than 0.2 eV have been presented 
by Rutherford & Vroom (1973), fitted by Linder et al. (1995) 
and then included in the early universe model calculations of 
Schleicher et al. (2008), but no direct calculations have been 
made available thus far. 

In this paper we therefore report an accurate analysis, from 
ab initio data, which deals with the above reaction. We present 
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new quantum calculations for process (2) over a broad range of 
collision energies: the corresponding rate coefficient has been 
obtained from the calculated reactive cross sections and is pre- 
sented in the following sections, where newly computed abun- 
dance of HeH + is also reported. Finally, we additionally com- 
pare it with the known behavior of LiH + cations under the same 
environmental conditions and can therefore draw some conclu- 
sions on their relative roles during the early universe chemical 
evolution. 



2. The quantum reactive calculations 

The elementary reactions indicated in either direction by pro- 
cess (2) have been studied both experimentally (Rutherford & 
Vroom 1973, Karpas et al. 1979, Tang et al. 2005) and the- 
oretically (Aquilanti et al. 2000, Palmieri et al. 2000, Xu et 
al. 2008, Ramachandran et al. 2009), because of its peculiar 
resonant features detected experimentally. Several potential en- 
ergy surfaces (PES's) have been obtained and employed, using 
quasi-classical trajectories, close coupling (time-independent 
and time-dependent) methods to study the inverse of the reac- 
tion (2). 

The polar cation of the present study, HeH + , has a large 
dipole moment (1.66 Debye) and has about 11 vibrational lev- 
els (Zygelman et al. 1998). For this reason it is therefore consid- 
ered of great importance in the early universe as possible coolant 
during the post-recombination era. As suggested by Dubrovich 
(1994), molecules with large dipole moments could interact with 
the cosmic microwave background radiation (CMB) inducing 
both spatial and spectral distorsions (Maoli et al. 1994). The 
interaction with the CMB is therefore crucial to determine the 
evolution and abundance of primordial molecules. The study of 
the HeH + chemical evolution, therefore, strongly depends on the 
accuracy of the main reaction rate coefficients which lead to its 
destruction or formation. Formation mechanisms for HeH + have 
already been extensively investigated (Zygelman et al. 1998) 
while we wish here to focus on its destruction mechanisms, and 
in particular on the direct exothermic (0.748 eV) reaction (2). In 
our calculations we employ one of the latest of the available, cal- 
culated PES (Ramachandran et al. 2009). Since it is also known 
that, under early universe conditions molecules are most likely 
to be in their roto-vibrational ground state (v = 0, j = 0), we 
have further decided to study reaction (2) starting from this level 
of the target cationic molecule. 

The exact quantum study of reactions involving an ionic sys- 
tem is usually rather difficult due to the presence of a long range 
"tail" in the interaction which therefore requires the use of a 
large basis set to generate numerically converged cross sections 
and rates. The use of an accurate computational method which is 
less costly in terms of resources is therefore mandatory. In a re- 
cent paper (Tacconi et al. 201 1) we have presented a novel com- 
putational implementation of the negative imaginary potential 
(NIP) scheme proposed earlier by Baer and coworkers (1990) 
and we shall use it for the present study. We start by solving 
the Schrodinger equation in a Body-Fixed (BF) frame of refer- 
ence, in the presence of an additional NIP. Hence we write the 
Hamiltonian of the triatomic system in the BF frame, within the 
coupled states (CS) (McGuire 1976) approximation 



2{i R 2fxR 2 2m r 2mr 2 

+V(r,R,-&) + V mp (r,R,$) (3) 



where R, r, and # are the Jacobi coordinates, fi and m are the 
reduced masses of the triatomic complex and diatomic molecule 
respectively. J and Q. are the total angular momentum and its 
projection on the relative axis while j is the diatomic rotational 
angular momentum. V(r,R,&) is the interaction potential and 
Vnip has the functional form given in Bovino et al. (2010b). The 
CS approximation provides the exact dynamics when the tar- 
get rotational angular momentum j is chosen to be equal to 
(MgGuire 1976). Thus, it is here equivalent to the full coupled- 
channel (CC) formulation. Because of the flux-absorbing effect 
coming from the NIP, the resulting S -matrix is now non-unitary 
and its default to unitarity yields here the (state-to-all) reaction 
probability, as discussed in Tacconi et al. (201 1). From the reac- 
tion probability one can in turn obtain the reactive cross section 

while the rate coefficients are computed by averaging the appro- 
priate reactive cross sections over a Boltzmann distribution of 
velocities for the incoming atom: 

a(T) = (8fc |f (f^ 2 f cr^. dU (E)exp(-E/k B T)EdE (5) 
(k B TY Jo 

where T is the gas temperature, k B is the Boltzmann constant 
and /i the reduced mass of the system. In the following section 
we shall present the effects from the accurate rates computed 
above on the final chemical evolution. 

3. The chemical network and the evolutionary 
modelling 

The evolution of the pregalactic gas is usually considered within 
the framework of a Friedmann cosmological model and the pri- 
mordial abundances of the main gas components are taken from 
the standard big bang nucleosynthesis results (Smith et al. 1993). 
The numerical values of the cosmological parameters used in 
the calculation are obtained from WMAP5 data (Komatsu et 
al. 2009) and are listed in Table 1 (for details, see Coppola et 
al. 2011). 

The abundance of HeH + is obtained by solving a set of dif- 
ferential coupled rate equations of the form 

—j£ = ffform - adest"l, (6) 

where af orm and a^est are formation and destruction rates for the 
reactant species ;' with number density n„ which in general de- 
pend on the number densities of the other species and on the radi- 
ation field. The rate coefficients involved in the helium chemistry 
network are the same as those reported in Galli and Palla (1998) 
(hereafter GP98), with the exception of the destruction reaction 
(2) which indeed now comes from the present quantum calcu- 
lations. The additional equations which govern the temperature 
and redshift evolutions have been reported in several earlier pa- 
pers (see e.g. GP98, Bovino et al. 2011) and therefore we shall 
not discuss them within the present context. 

4. Results and comparisons 

As already outlined in the Introduction, our understanding of 
the early universe evolution strongly depends on our knowl- 
edge of the elementary processes involving atomic and molec- 
ular species taken to be present at that time. At the molecular 
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Table 1. Cosmological model parameters. Ho is the Hubble con- 
stant, Q m , £2 r , Q A , Q K are density parameters, and fi are the ini- 
tial fractional abundances of the main species. To is the present- 
day CMB temperature (see Coppola et al. 201 1 for details). 



Parameter 


Numerical value 


H 


70.5 km s-' Mpc- [ 


T 


2.725 K 


Zeq 


3141 


fidm 


0.228 


fib 


0.0456 


n m 


fidm + fib 


fi r 


fim/(l + Zcq) 


fi A 


0.726 


C2k 


1 — £2 r — f2 m — 


/h 


0.924 


/He 


0.076 


/d 


2.38 x 10~ 5 


fu 


4.04 x 10- 10 



level, therefore the observables that one needs to obtain are the 
various cross sections related to their mechanisms of occurrence. 
In the present work we have accurately calculated the cross sec- 
tions for HeH + in reaction with H( 1 5) using quantum dy- 
namics and ab initio data which provide the interaction forces. 

In the following we shall report and discuss our results and 
further compare the new HeH + abundances produced by the new 
quantum reaction data with those obtained for another molecular 
ion present in the chemical network, LiH + (Bovino et al. 201 1). 

4.1. The reaction cross sections 

A set of three hundreds coupled equations which describe re- 
action (2) have been solved using our recent, in-house code 
(Tacconi et al. 2011). A direct product of a Colbert-Miller 
(Colbert & Miller 1992) discrete variable representation (DVR) 
of 100 points ranging from 0.5 «o to 12 do for the reagent ion 
vibrational coordinate and a rotational basis of 40 spherical har- 
monics has been employed to describe the reactive complex. The 
solution matrix has been propagated up to 50 do and the con- 
vergence was checked as a function of the propagator step-size, 
yielding a final convergence of cross section values within about 
1%. 

The calculations were carried out over a wide range of en- 
ergies, from 10~ 4 to 1.0 eV, and for total angular momentum 
values (J) ranging from 10 to 90 for the highest energy. The 
NIP parameters were tested following the Baer criteria (Baer et 
al. 1990): the NIP stability has been reached for r m ; n = 2.75«o, 
'"max = 5.75«o and the NIP order n - 2 (see Tacconi et al. 201 1). 

Our computed, /-converged integral cross sections are pre- 
sented in Figure 1 . It is clear from the Figure the non-Langevin 
behavior of the cross sections for this ionic reaction since, in- 
stead of becoming independent of E at lower energies, they go 
through a marked maximum around 10~ 2 eV. In the same fig- 
ure we report the experimental (thick line, from Rutherford et 
al. 1973) and the extrapolated (dotted line, from Linder et al. 
1995) cross sections which reach down to collision energy val- 
ues around 10 -2 eV. The agreement between our theoretical re- 
sults and the experimental findings, both measured and extrapo- 
lated, is reassuringly good and confirms the physical reliability 
of the present NIP method for handling ionic reactions. No mea- 
sured data are available at the lower collision energies as crossed 
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theoretical (present calculations) 

experimental {Rutherford et al.) 

extrapolated (Linder et al.) 



Energy (eV) 



Fig. 1. Computed integral cross sections as a function of en- 
ergy. The thick line represents the fitted experimental data from 
Rutherford & Vroom (1973). The dotted line is obtained from 
the formula given by Linder et al. (1995) by extrapolating the 
Rutherford's data. 
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Fig. 2. Computed production/destruction rates of HeH + as func- 
tion of redshift. Solid and dashed curves represent production 
and destruction processes, respectively. 



beam experiments with ionic partners become increasingly dif- 
ficult in that range. It is, however, worth noting here that the 
results of Figure 1 are very different from the findings for LiH + 
(Bovino et al. 2010a) in the same range of energies. Such differ- 
ences are related to different microscopic mechanisms of the re- 
action pathways between the two systems, an aspect of the prob- 
lem which we are planning to analyse more in detail in further 
work. Qualitatively, the two systems, although both undergoing 
no barrier, exothermic reactions, have very different exothermic - 
ity values that affect product separations at low energy. 
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Table 2. Computed rate coefficient for the reaction (2) as a func- 
tion of T 



1(K) 


a (cm 3 .? ') 


i 
1 


A A 1 1 r\— 12 
4.41 X 10 


c 

5 


2.1 1 X 10 


10 


c Ai w i r\— 1 1 

5.4/ X 10 


20 


1.33 x 10 






50 


3.43 x 10-'° 


70 


4.39 x 10" 10 


80 


4.76 x 10-'° 


100 


5.33 x 10-'° 


200 


6.76 x 10" 10 


300 


7.39 xlO- 10 9.10 + 2.5 x 10" 10 (*) 


500 


8.08 x 10" 10 


800 


8.70 x 10-'° 



2500 400 100 40 10 1 



a 





LiH + + H — > Li + + H2 




calculated 

: 


HeH + + H — * ri, + + He : 




- — calculated 




— extrapolated {Schleicher et al.) 




• — • experimental (Karpas et al.) 





I I I I I I I 

1000 100 10 1 



Temperature (K) 



(*) experimental (Karpas et al. 1979) 

4.2. The depletion rates 

The new rate coefficients are reported in Table 2 and further fit- 
ted with the following formula: 

a(T) = 4.3489 x 1(r 10 r (m037 V- 31 - 53%/r) (7) 

valid for T < 1000 K. The experimental value obtained by 
Karpas et al. (1979) is also reported for comparison in the same 
Table. It is again very reassuring to see that the only existing 
experimental datum at 300 K is indeed very close to our compu- 
tational value at the same temperature. 

In order to more clearly evaluate the significance of the 
newly computed cross sections we need to view the process (2) 
within the broader context of the variety of elementary events 
involved with this molecular cation. To that end, Figure 2 shows 
the main processes leading to the formation and the destruction 
of HeH + . As clearly shown by the figure, the radiative associ- 
ation process (pi) dominates the formation of HeH + over the 
whole range of redshifts, while the association involving the He 
ion (p2) and the one involving the HJ molecule (P3), have negli- 
gible contributions at the lower redshift values, the most impor- 
tant range for the present study. 

From the same figure one further realizes that the photodis- 
sociation mechanism (d\) is the dominant destruction path of the 
HeH + for z > 300 while, on the other hand, the chemical path 
which involves the reaction (2), the one which destroys HeH + 
(c/2) in collision with H, is seen to be the most important pro- 
cess for z < 300. A contribution also comes from the electron- 
assisted dissociation reaction (J3) that neutralizes the charged 
partners and dominates over (di) for z < 10. 

4.3. A comparison with LiH + 

We do know that during the Dark Ages of the universe, the exist- 
ing molecules interacted with CMB photons and with the present 
atoms. What we do not know is if these species existed with large 
enough abundances to produce now detectable signals. Because 
of its large dipole moment (1.66 Debye), and the large abun- 
dances of its atomic components (H and He), HeH + , is certainly 
one of the prime candidates amenable to observations (Persson 
et al. 2010). This is even more important as our recent ab-initio 
results (Bovino et al, 2009, 2010a, 2010b, 2011) showed that 
lithium-containing molecules such as LiH and LiH + , in spite of 



Fig. 3. Depletion rate coefficients as a function of the tempera- 
ture. In the upper scale the redshift is also reported. The experi- 
mental rate coefficients of Karpas et al. (1979) and those extrap- 
olated by Schleicher et al. (2008) are also presented. The lithium 
reactions rate is from Bovino et al. (2010a). 



their significant dipole moments that would favor transitions (5.9 
and 0.6 Debye), have very low abundances over the whole range 
of redshifts studied. 

In Figure 3 we report the calculated rate coefficients for the 
depletion reaction (2) over a wide range of temperatures (red- 
shifts) together with the earlier calculations for LiH + depletion 
(Bovino et al. 2010a). As it is clear from that Figure, the HeH + 
+ H reaction turns out to be very efficient at the redsfhifts above 
100, with rate values comparable to those of the LiH + + H de- 
struction reaction given in the same figure (dot-dashed line). 
However, while the latter remains essentially constant over the 
whole range of examined redshifts, the efficiency of reaction (2) 
quickly decreases for z < 100, coming down by about two orders 
of magnitude by the time unitary redshift values are reached. 
This behavior therefore favors the HeH + as a more likely sur- 
vival candidate vis-a-vis the LiH + case, thereby providing a good 
molecular coolant acting during the primordial star formation. 
In the same figure we also report the experimental rate at 300 K 
given by Karpas et al. (1979); the rates included in the GP98 
model are taken from the experimental data of Rutherford & 
Vroom (1973). As seen earlier for the cross sections, the ex- 
perimental rates are in good agreement with our quantum results 
(see also Table 2). 

In Figure 4 we report the calculated abundances for HeH + 
obtained using our new quantum rates discussed above: they are 
compared with the old ones obtained by GP98 and also with 
our newest results for LiH + (Bovino et al., 2011). The scenario 
emerging from the above calculations indicates that for redshift 
around z - 20, the molecular HeH + abundance has now in- 
creased by more than one order of magnitude with respect to 
the earlier calculations. It thus remains a more abundant ionic 
species LiH + , hence becoming an interesting molecular candi- 
date for actual detection. As shown by Maoli et al. (1994) (see 
also Schleicher et al. 2008), the relevant quantity for determin- 
ing the imprint of molecules on the CMB is the optical depth due 
to line absorption. In the case of HeH + , the optical depth has a 
broad peak around a frequency of ~ 50 GHz, due to redshifted 
rotational transitions with j = 4-6 of the ground vibrational 
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Fig. 4. Computed relative abundances of HeH + in the post- 
recombination era as function of redshift z: present results (solid 
curves), results of GP98 (dashed curves). The LiH + results are 
from Bovino et al. (201 1). 

level. With the new rate coefficient derived in this work, the max- 
imum value of the optical depth, ~ 10~ 7 , is about one order of 
magnitude larger than the value computed earlier by Schleicher 
et al. (2008), an encouraging result for observational perspec- 
tives. The new reaction rate has also an effect on the abundance 
of , that is reduced by a factor of ~ 5 in the range 10 < z < 50 
with respect to earlier results. However this reduction has no sig- 
nificant effect on the abundance of H 2 because the Hj channel 
of production of H2 is no longer active at these redshifts. 

5. Conclusions 

We have carried out a new set of ab initio, quantum calculations 
for the rate coefficient of the reaction (2), which therefore per- 
mit us a revision of the abundance of HeH + molecule within the 
expected conditions operating during the early universe. One of 
the main results from the present reactive calculations is that the 
chemical destruction rate of this polar cation exhibits a tempera- 
ture dependence that turns out to be very different from previous 
data, both experimental and theoretical, especially for tempera- 
tures smaller than 100 K. 

These new results obtained for the HeH + abundance now in- 
dicate that this molecule should be much more likely to have 
survived at low redshift, since its fractional abundance goes up 
by more than one order of magnitude from previous estimates, 
becoming of the orderof 10 . This new findings should be con- 
sidered within the broader context whereby the fractional abun- 
dance of LiH + , recently calculated by Bovino et al. (201 1), turns 
out to be significant mostly in the low-redshift region of z < 30. 
The comparison between the specific fractional abundances of 
the two species therefore suggests rather compellingly that the 



HeH + ionic molecule, if more abundant than the LiH + , would 
become more likely to be observed experimentally, even more 
so than what was suggested by previous works (GP98, Persson 
etal. 2010). 

It is also worth to mention that, due to its larger dipole mo- 
ment and to its greater number of roto-vibrational levels, when 
compared with LiH + , the molecular cation HeH + turns out to be 
the best candidate to date for experimental observations in dif- 
ferent astrophysical environments, a result now supported by the 
present ensemble of accurate computational studies from quan- 
tum, ab initio methods. 
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